Poynting's theorem and energy conservation 
in the propagation of light in bounded media 
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Starting from the Maxwell-Lorentz equations, Poynting's theorem is reconsidered. The energy 
flux vector is introduced as S e = (E x B)//io instead of E x H, because only by this choice the 
energy dissipation can be related to the balance of the kinetic energy of the matter subsystem. 
Conservation of the total energy as the sum of kinetic and electromagnetic energy follows. In 
our discussion, media and their microscopic nature are represented exactly by their susceptibility 
functions, which do not necessarily have to be known. On this footing, it can be shown that energy 
conservation in the propagation of light through bounded media is ensured by Maxwell's boundary 
conditions alone, even for some frequently used approximations. This is demonstrated for approaches 
using additional boundary conditions and the dielectric approximation in detail, the latter of which 
suspected to violate energy conservation for decades. 

PACS numbers: 03.50.De, 78.20.-e, 77.22.Ch, 41.20.Jb, 71.35.Cc 
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INTRODUCTION 

The effects of a medium on radiation (e.g., in an ex- 
periment on the transmission and reflection of light in a 
dielectric medium) are determined exactly by the suscep- 
tibility function x(r, r', t, t', E), which enters Maxwell's 
equations through the material equations (also known as 
constitutive relations). Due to its complexity, it is only 
in rare cases that \ can be expressed analytically or valu- 
able relations can be drawn from it. For a long time, the 
more complex properties of \ have not been of much in- 
terest, and textbooks still introduce the susceptibility as 
a constant scalar or as a constant tensor, to account at 
least for anisotropy (e.g., [H, Q). Non-linearity (i.e., de- 
pendence on E) is an own interesting field of research and 
technology but out of scope here and will be neglected. 

But strictly spoken, the presence of a medium bound- 
ary alone enforces a non-local and spatially inhomoge- 
neous x, if spatial dependence in the dispersion rela- 
tion (spatial dispersion) is not explicitly neglected. In 
crystal and semiconductor optics, some approaches to 
this problem resort to using the susceptibility of an in- 
finite medium for the description of a bounded real one. 
Then, multiple polariton modes are predicted to propa- 
gate concurrently in the medium, and the precondition 
of a continuous transition of the electromagnetic fields at 
the medium boundaries (Maxwell's boundary conditions, 
MBC's) is no longer sufficient to determine the ampli- 
tudes and phases of these modes. As a solution, dif- 
ferent phenomenological additional boundary conditions 
(ABC's) were proposed d Q, the first of which in 1957. 
In 1973, a related approach was published 5|, suggesting 
to limit the bulk susceptibility to the bounded medium 
by step functions (dielectric approximation, DA). These 
methods were applied in a vast number of publications 
in the past decades, and they still are used today and 
sucessfully explain experimental data 



However, it is known that some of these approxima- 
tions can conflict with the fundamental law of energy 
conservation 0, Q. Thus, scientists were prompted to 
take greatest care in their argumentation or even to re- 
ject these approximations despite of their successes (see, 
e.g., @,i). 

We will show that, in the case of the DA and other 
important approximations, energy conservation is not a 
question of the approximation but rather is ensured by 
the MBC's alone. The DA is addressed in detail: The 
compatibility of the approximation with energy conser- 
vation is proven, and problems in the respective work are 
pointed out. For all this, we use an approach that one 
of us recently developed [lj| , which provides analytically 
exact relations between radiation and matter in the elu- 
sive case of spatial dispersion with the aid of Keldysh's 
photon Green's functions. 

Throughout the discussion of bounded media, we re- 
gard arbitrary steadily excited media. They are repre- 
sented by their susceptibility functions, which do not nec- 
essarily have to be known. Besides the latter, the corner- 
stone of our discussion is Poynting's theorem (PT), the 
equivalent of the energy conservation law in electrody- 
namics. It balances the time derivative of electromag- 
netic and mechanical energy U e>m and sources of respec- 
tive energy flux S e , m : 



U„ 



div (S e 







(1) 



In Poynting's original work llj and many established 
textbooks [l|, the theorem is given — without re- 

garding spatial dispersion — as 

ED + HB + div S e = -j ext E, S e = div (E x H). (2) 

An expression for the electromagnetic energydensity U e 
can only be given with further assumptions [2 1| . This is 
why eq. (|2|) cannot be generally interpreted as an energy 
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continuity equation and, consequently, div (E x H) not 
as an energy flux vector. 

In order to lay solid grounds for our considerations, 
we will at first give a thorough derivation of the theorem 
((T|) in conjunction with spatial dispersion and mechanical 
energy, giving hints on why several articles, such as [l3| 
or those referenced therein, doubt the validity and form 
of eq. (J2j) in more complex situations. 

A note on terminology: Historically, the term Poynt- 
ing vector denoted the electromagnetic energy flux S e 
only. With regard to eq. |T]), we here need to distinguish 
clearly between the electromagnetic and the mechanical 
Poynting vector. 



A CLEAN-ROOM DERIVATION OF 
POYNTING'S THEOREM 

There are many doubts expressed on under which con- 
ditions (static magnetic fields, stationary situations, dis- 
persive media) PT applies and on how Poynting's energy 



flux vector has to be interpreted (compare [13. 1 14 . Il5l . 
[lij). Here, we do not see any reason for these. Instead, 
we will give a brief derivation of PT, while carefully ac- 
counting for any approximations made or conditions im- 
plied. This derivation is done in classical physics but 
can easily be extended to quantum physics by applying 
the well-known recipes (i.e., operators, symmetrization, 
commutation relations). 

Surprisingly, we do not need any assumptions about 
the matter. Avoiding constitutive relations and derived 
fields, we rather rely on the microscopic Maxwell-Lorentz 
equations [I7j . which are the most fundamental and uni- 
versal: 

rot E = B, rot B = + -4-E, (3a) 
j{r,t) =^2e Q r a 6{r - r a (t)), (3b) 

a 

with j as a current of charged particles a. 

With the vector identity div (X x Y) = Y rot X — 
XrotY, we immediately have 

e EE + — BB + —div (E x B) = -jE (4) 

Mo Mo 

and may identify the electromagnetic field energy density 
and the electromagnetic Poynting vector as 



i( £o E 2 + ^B 2 ), 



Mo 



1 



S e = — (E x B) 

Mo 



(5) 
(6) 



Note that in what follows, we will use these new defini- 
tions of the two quantities instead of the more conventinal 
given in eq. [[2]). if not otherwise stated. 



Now, a microscopic mechanical expression for jE shall 
be developed. We use the density of mechanical (kinetic) 
energy of the particles a, U m , as well as their kinetic 
energy flux, S m , according to 



CUr,t) = £^r^(r-r a (i)), (7) 

a 

S ro (r,t) = £^r2*a*(r-r a (i)). (8) 



The force mf appearing in the time derivative of U m is 
just the Lorentz force F = e(E + r x B), yielding 



d 

— U m (r,t) +divS m (r,i) 



:j(r,t)E(r, t). 



(9) 



This can finally be combined with eq. @ to Poynting's 
theorem as shown in eq. |T]). No restrictions have been 
made to the behaviour of the fields or particles or the 
nature of the latter, except for non-relativistic velocities. 
This approach will turn out to be inevitable for estab- 
lishing a relation between electromagnetic and mechani- 
cal energy. Mechanical energy here is purely kinetic, any 
potential energy is attributed to the fields. 

The continuity of the energy flux at medium bound- 
aries is an aspect to which much attention is paid in 
literature. The electromagnetic Poynting vector S e de- 
rived here is not necessarily continuous, but this does 
not account for physical meaningfulness. It is the energy 
continuity eq. Jl} that holds. The continuity of one the 
quantities involved depends on that of the others, which 
in turn is subject to the system and model considered. 



MATERIAL EQUATIONS AND 
SUSCEPTIBILITY 



The macroscopic Maxwell equations originally were in- 
troduced to obtain a global, effective description of elec- 
trodynamics in the presence of a medium. One may con- 
sider their form a well-founded historic convention, which 
can be reproduced simply by splitting j = jj n( j + joxt and 
the ansatz 



jind = P + rot M, 



(10) 



which is also suggested [18[ by the well-known multi- 
pole expansion of the electromagnetic properties of the 
medium and the respective material equations 



D(r,i) 


= e E(r,t)+P(r,i) 


(11a) 


H(r,t) 


= — B(r,t)-M(r,i) 


(lib) 


Mo 


P(r,i) 


= e Q J dVdt' Xcicc(r,r',i,i')E(r',i') 


(11c) 


M(r, t) 


= J_ /"dVdt' Xmagn(r,r',M')B(r',t')- 
Mo J 


(lid) 
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Note the clear analogy in this notation and the ob- 
vious distinction of fields E,B and derived fields D,H. 
Usually, textbooks [J, 0, [l8| focus more on the analogy 



B — /!()ArH, D — £(]£ r E, 



(12) 



i.e., with constant tensorial relations (|llc|) . (|lld[) , and 
consequently treat E and H on the same level. But the 
physically relevant fields are actually E and B, because 
they determine, e.g., the Lorentz force, the energy den- 
sity, and the electromagnetic energy flux. 

Neither the multipole expansion discussion in [12j nor 
the tensorial descriptions in [l], Q end up with non-local 
or inhomogencous material equations and thus cannot be 
used with spatial dispersion and bounded media. Hence, 
a more general dependency of the derived fields on E 
and B, as shown in eqs. HJ), needs to be used. Here, 
the tensorial susceptibilities x c icc, Xmagn represent ex- 
actly the linear properties of an arbitrarily complex sta- 
tionary medium. However, PT can still be derived in 
the same manner without any restrictions to the electro- 
magnetic fields, properties of the medium or mathemat- 
ical properties of Xcicc or Xmagn, except for the temporal 
homogeneity of the latter. Also, eq. ([2]) can still be ob- 
tained, but then j does no longer fulfil eq. (|3bj) and thus 
cannot be coupled with the mechanical energy continuity 
relation @, and the field energy density cannot be given. 

In optics, the approximation Xmagn = is widely ac- 
cepted. Then, S e becomes E x H and j in d = P. From 
a higher level of abstraction one can also postulate that 
jind = P- This avoids pinpointing magnetization and 
other difficulties and results in the following material 
equations 19]: 



D(r,t)=e„E(r,t)+P(r,t) 

P(r, t) = e J dVdi' X (r, r', t, t')E(r', t') 



H(r,i) = — B(r,i) 



(13a) 
(13b) 

(13c) 



Now, the H field is merely a rescaled B field and any 
electromagnetic effects are contained in P. 

Much of the controversy around the validity of Poynt- 
ing's theorem may be due to inconsistent formulation of 
the material equations or using eq. © despite of magne- 
tization. In 13j, the author correctly finds that eq. ([2]) 
together with the common approximation (|12[1 for the 
material equations does not apply for wide categories of 
interactions, such as spatial dispersion. He also bases his 
further considerations on the Maxwell-Lorentz equations 
and equally finds eq. ([4]) but then develops an elaborate 
expression for jE within the (local) multipole expansion 
approximation. The result, of course, is an expression for 
our S m in eq. UJ (compare [H, eq. 21]). This is an hon- 
orable effort, even though his results can only hold true 
in the limits of this approximation. Unfortunately, the 



author does not identify (E x B)//xq as a valid electro- 
magnetic Poynting vector and eq. ([1]) as a more general 
energy conservation law. 



THE WORK OF BISHOP/MARADUDIN AND 
ITS CONCEPTUAL PROBLEMS 

In [tJ , Bishop and Maradudin analyze the energy flow 
in a semi-infinite spatially dispersive dielectric (i.e., half- 
space geometry) on which light is incident from the vac- 
uum. The energy flux in the vacuum is given by the elec- 
tromagnetic Poynting vector S e only. In the medium, 
the mechanical Poynting vector S m appears additionally. 
They state that the DA fails to conserve energy by show- 
ing that S m does not evolve continuously from zero at the 
surface, as if there was an energy source. For this pur- 
pose, they develop a model of the transport of mechanical 
energy in a special system, a diatomic cubic crystal, on 
the basis of fixed, undamped and uncoupled harmonic os- 
cillators only, then derive PT from Maxwell's equations 
and adapt it for their model. In 0, eq. 4.7] they postulate 
that S m be represented by the oscillations; the spatial re- 
striction to the medium is forced by an appropriate step 
function. 

This postulate does not necessarily fulfill the continuity 
condition. S m is continuous only if the oscillator ampli- 
tude vanishes at the surface. One of the possibilities to 
ensure this is Pekar's ABC [3J, as the authors confirm in 
their conclusion. It is no surprise that the step function 
leads to a delta-shaped term [7J, eq. 4.10] when conduct- 
ing the divergence operation in eq. |T]). Of course, such 
a source term is not compatible with energy conserva- 
tion. Consequently, the authors generalize their model 
and derive additional conditions for its free parameters. 

This kind of approach is highly doubtful because any 
predictions might have to be attributed to the model and 
its deficiencies. Here, this is clearly the case: The me- 
chanical Poynting vector as postulated does neither con- 
tain all mechanical energy transport in the medium (e.g., 
via phonons, center-of-mass motion of excitons etc.) nor 
does it fulfill itself the condition to be analyzed. 



COMPATIBILITY OF DA AND PT 

We rather suggest an approach free of models. It in- 
cludes but is not limited to the case of DA; it does not 
make any assumptions on the microscopic nature of the 
medium but instead is based on the susceptibility \ as 
an exact but unknown representation of the latter. 

We regard a stationary isotropic medium in the slab 
geometry. Light is incident perpendicular to the medium 
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surface. Then Poynting's theorem reduces to: 
dS e (x,t) dS m (x,t) 



dx 



dx 



j(x,t)E(x,t) =: -W(x,t) 
(14) 



By expressing E and B through the vector potential 
A(x, t), j through dP/dt and Fourier transforming t — > 
01, we obtain 



W{x,w)=ie J ^-J\lo-J) 



dx'\(x, x' ,uj')A(x' , ui')A(x, 01-01'). 



(15) 



There is no need for taking a time average. As x = 
in vaccuum, it may be replaced by x( x -> x ' \ w ) = ©("j — 
|#|)G(4 — \x'\)x(x,x' ,oj), where L is the length of the 
slab. Of course, this is also valid for the special case of 
DA, where x just simplifies to 0(^ — |x|)0(4 — \x'\)x(x — 
x',oi). Evaluating the step functions, we have: 

W(x,Lo)=is a I ^' 2 (cj-oi')e(^-\x\)A(x,o;-u') 



dx' x(x, x , ui)A(x , to') 



,L 



=:Q(--\x\)W(x,o>) 

On the left surface, the step function reduces to 0(|n + 
a;). We now obtain a new, model- free relation for S m by 
integrating (fl"4")) : 

S m (x,u) = Q(^ + x)J dx'W(x',oi) (16) 

As can clearly be seen, S m is zero at the surface and grows 
continuously for x > — -j, because W is a finite value. 
It is therefore continuous at the surface. There is no 
evidence for any energy source or violation of Poynting's 
theorem. 



ENERGY CONSERVATION WITH BULK 
SUSCEPTIBILITY AND ABC'S 

Notwithstanding the first successful microscopic non- 
local computation of the polarisation [20| , which required 
supercomputers, the pragmatic assumption of a bulk sus- 
ceptibility xili^) is still widely used to analyze light 
propagation through bounded media. In this case, the 
dispersion relation 



q(o>) 2 = -5-(l+x(9.w)) 



(17) 



can have various solutions qi(oj), giving rise for different 
polariton modes. Consequently, the electromagnetic field 



inside is a superposition of these forward and backward 
propagating modes (we stick with slab geometry): 



E(x,u)=Y;[f^) eiqi{U)X + h ^> 



-\qi {uj)x 



The extinction theorem, 
a\(oj) = 1 



(18) 



(19) 



where a\ is the absorption and r, t the coefficients of re- 
flected or transmitted light, resp., and the MBC's, 



e -iqo 2 -)_ T-e' 90 T 

g-igo-j- _ re igo-| 

te [qo ^ 
te i?0 * 



E(/« 



91 2 +bie 1Qi 2 

i 



(20a) 
(20b) 
(20c) 
(20d) 



where qg = lu/c, form a set which could be solved if only 
a single mode was present but in fact is under-determined 
due to the polariton modes. The postulation of ABC's 
help solving this problem, but some of them were known 
to conflict with energy conservation. (For a concise sum- 
mary on ABC's and the controversy around them, see § 
and refs. therein). Most of these ABC's, as well as the 
DA, can be reproduced by generalized ABC's (l9j |. 

Using the RHS of the integral representation of 
eq. (|14p . an independent expression for the absorption 
can be derived: 



a 2 (oi) 



1 01 

2i~ 



dxdx' A*(x, uf)x(x, x', o>)A(x',ui), (21) 



with x(x,x',uj) = x(x,x',w) - x*(x',x,oi) (22) 

It follows by applying a monochromatic wave and some 
straightforward calculation. Then, static terms (oc S(oi)) 
and such oc 5(oj ± 2wo) appear. The former can be re- 
garded seperately, elegantly circumventing cycle averag- 
ing. The important relations (fT9| and (f2Tj) are presented 
in [l 01 ] in a quantum mechanical generalization valid for 
non-locally dispersive media and oblique incidence. They 
take almost the same form as shown here. 

The polarisation here follows from eqs. (flU)) as 



P(x,oi) = £xiM hiM&M* + biMe-*™* 



It can be used to calculate 0-2(01) in a variant of eq. ^1 
based on field strength and polarisation function only, 



02(01) = — Im 
c 



dxE*(x,oi)P(x,o>). 



(23) 
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Now, constructing eq. (j 19[) as 

ai (o,) - i( KEOH) + d20bl)| 2 - |(|20iD - dSoEl)! 2 

— jdSQcJ) + ([20d]>| 2 + KHOSD — <[2Qd!)| 2 ), 

the last term equaling zero, it can be shown that a\ =02- 
As no ABC's have been used up to now, this means that 
energy conservation in such a system is ensured by the 
MBC's alone. 



of Pekar's and a dead layer, see [9|). 

All these results save a plethora of recent and older 
work the shortcoming of possibly violating fundamental 
principles, and the use of their results in research is finally 
justifiable, at least in this respect. 

The authors would like to thank the Deutsche 
Forschungsgemeinschaft for support through Sonder- 
forschungsbereich 652. 



CONCLUSION 

We hopefully could present the concept of Poynting's 
theorem in conjunction with a mechanical energy flux 
vector in an illustrative way, concretizing and extending 
textbooks in this sense, offer a best-practice approach, 
and show that doubts on its validity in certain conditions, 
such as spatial dispersion, arc unsubstantiated. By the 
given approach, a clear relation between electromagnetic 
and kinetic energy in a system can be established, and 
electromagnetic energy density and flux are determined 
by the fundamental physical fields E and B alone. 

We then have proven analytically that energy conser- 
vation in the propagation of light through bounded media 
is fulfilled for arbitrary susceptibilities, and thus also for 
the dielectric approximation. 

Often, the use of models and approximations for the 
microscopic nature of a medium is misleading, and any 
predictions have to be taken with a grain of salt. This 
was pointed out as a problem in Q and in the common 
perception of Poynting's theorem in general. Great care 
must be taken in the application of the material equa- 
tions, and it seems advisable to always start from (fTTj) 
and the Maxwell-Lorentz equations j3|), especially when 
spatial dispersion is dealt with. If a model for the me- 
chanical Poynting vector is applied, it must be consistent 
with the susceptibility. Namely, it can be constructed via 
eq. (fT6]). 

Furthermore, we were able to prove analytically that 
energy conservation in approximations employing bulk 
susceptibility is completely assured by Maxwell's bound- 
ary conditions alone, even if multiple polariton modes are 
predicted. Contrary to what was thought before, addi- 
tional boundary conditions do not play any role in this re- 
spect, whether positive or negative. Consequently, from 
the energy conservation point of view, it is possible to 
freely choose the ABC's that match microscopic calcula- 
tion results best (which would probably be a combination 
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